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We analyze the fluctuating dynamics at the golden-mean transition to chaos in the critical circle 
map and find that trajectories within the critical attractor consist of infinite sets of power laws 
mixed together. We elucidate this structure assisted by known renormalization group (RG) results. 
Next we proceed to weigh the new findings against TsaUis' entropic and Mori's thermodynamic 
theoretical schemes and observe behavior to a large extent richer than previously reported. We find 
that the sensitivity to initial conditions has the form of families of intertwined g-exponentials, 
of which we determine the g- indexes and the generalized Lyapunov coefficient spectra \q. Further, 
the dynamics within the critical attractor is found to consist of not one but a collection of Mori's 
g-phase transitions with a hierarchical structure. The value of Mori's 'thermodynamic field' variable 
q at each transition corresponds to the same special value for the entropic index q. We discuss the 
relationship between the two formalisms and indicate the usefulness of the methods involved to 
determine the universal trajectory scaling function a and/or the ocurrence and characterization of 
dynamical phase transitions. 
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I. INTRODUCTION 

The unfamiliar dynamics at critical attractors of even 
the standard one dimensional nonlinear iterated maps 
0-0 is of current interest to statistical physicists; it 
provides insights into properties of systems at which the 
key supports of the canonical theory, phase space mix- 
ing and ergodicity, break down. At these attractors the 
indicators of chaos withdraw, such as the fast rate of sep- 
aration of initially close-by trajectories. There is a small 
number of known routes to chaos and it is of interest to 
know if at their thresholds the dynamics displays simi- 
lar features and if so to what extent. If present, these 
similarities would apply above the known universality of 
each type of critical attractor. 

Critical attractors of nonlinear low-dimensional maps 
have always drawn much interest because they demon- 
strate the ways by which nonlinear systems bridge peri- 
odic and chaotic motion The geometrical prop- 
erties of these attractors have been known for a long 
time but their dynamical properties are even now 
perplexing and still need to be full y u nderstood. Re- 
newed interest in this subject |^ - |l5l | has focused on 
the applicability of g-statistics, the generalization of the 
Boltzmann-Gibbs (BG) statistical mechanics proposed 
by TsaUis |0, to these nonlinear systems. Impor- 
tant examples of critical attractors are the onset of chaos 
via period doubling, intermittence and quasiperiodicity, 
the three universal routes to chaos in dissipative maps 
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displayed by the prototypical logistic and circle maps [ij 
- ^3]. Perhaps the most interesting types of critical at- 
tractors are critical strange nonchaotic attractors as 
in the period doubling and quasiperiodic onset of chaos. 
These are geometrically complicated multifractal attrac- 
tors with a vanishing Lyapunov coefficient Ai and a sen- 
sitivity to initial conditions that does not converge to 
any single-valued function but instead displays a fluctu- 
ating pattern that grows as a power law in time t - 

M 

Trajectories within a critical attractor show self-similar 
temporal structures, they preserve memory of their pre- 
vious locations and do not have the mixing property of 
truly chaotic trajectories A special version of the 

'thermodynamic formalism' was adapted long ago j2nl |- 
[23| to study the dynamics at critical attractors and 
quantitative results were obtained that provided a first 
understanding, particularly about the envelope of the 
fluctuating and the ocurrence of a so-called 'g-phase' 
dynamical phase transition - As shown below 

we uncover a complete and rich dynamical structure for 

and clarify its relationship to families of g-phase transi- 
tions. We explain how these transitions relate to distinct 
regions of the multifractal attractor, and their manifes- 
tation in the universal jump discontinuities of the trajec- 
tory scaling function a |24j - . We clarify also the rela- 
tionship of the thermodynamic approach with g-statistics 
with regards to the dynamical properties studied. 

Recently, rigorous results have been developed 
that support the validity of g-statistics (outlined in the 
next paragraph) for the critical attractors associated with 
the intermittency and period-doubling routes to chaos, 
i.e. the tangent bifurcation and the accumulation point 
of the pitchfork bifurcations, but studies of the same type 
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have not been carried out for the quasiperiodic route to 
chaos. Here we present results on this route by consid- 
ering the proverbialgolden-mean onset of chaos in the 
critical circle map 0|, Q and show that there are el- 
ements of universality in their dynamical properties as 
they appear in close correspondence with those recently 
analyzed for the Feigenbaum attractor [T5j |. 

The central point of the (/-statistics with regards to 
the description of the dynamics of critical attractors is 
a sensitivity to initial conditions related to the q- 
exponential functional form, i.e. the 'g-deformed' expo- 
nential function exp,(a;) = [1 - (q - l)a;]-i/(9-i) [|. 
From such ^( one or more spectra of g-generalized Lya- 
punov coefficients A, can be determined 0. The 
are dependent on the initial position and each spectrum 
can be examined by varying this position. The Xq satisfy 
an identity Xq — Kq where Kq is an entropy production 
rate based on the Tsallis entropy Sq, defined in terms of 
the q-logarithmic function \n„ y = (?/^~^— 1)/(1 — g), the 
inverse of exp^(a;) 3, 0, [Tsj . 

Here we find that at the golden-mean onset of chaos 
takes the form of a family of interwoven g-exponentials. 
The g-indexes appear in conjugate pairs, q and Q = 2 — q, 
as these correspond to switching starting and finishing 
trajectory positions. We show that q and Q are related 
to the occurrence of pairs of dynamical q-phase transi- 
tions that connect qualitatively different regions of the 
attractor - These transitions are identified as 

the source of the special values for the entropic index 
q, and these values are determined by the universality 
class parameters to which the attractor belongs. In our 
present case the parameters are the golden-mean winding 
number Wgm = {V5 — l)/2 ~ 0.618034 and the universal 
constants (in the sense of Feigenbaum's trajectory scaling 
function a |2^ - pG*]) that measure the power-law clus- 
tering of iterate positions. The most prominent of the 
universal constants is agm — —1.288575, that appears in 
the description of the most rarefied region and (as a^m) 
of the most crowded region of the multifractal attractor. 

We are interested here in determining the detailed de- 
pendence of the aforementioned dynamical structure for 
the golden-mean chaos threshold on both the initial po- 
sition and the observation time t as this dependence is 
preserved by the infinite memory of trajectories. To this 
purpose we recall in the following Section 2 the famil- 
iar features of the circle map, and the properties of the 
renormalization group (RG) fixed-point map for a cubic 
inflection point with zero slope. Then in Section 3 we 
describe the organization of the superstable trajectory 
at the chaos threshold into families of positions that lie 
along power laws in time. Assisted by these properties, 
in Section 4 we determine the sensitivity for trajecto- 
ries that have their starting and finishing positions at the 
most crowded and most sparse regions of the multifrac- 
tal attractor. In Section 5 we explain how the properties 
for and its Xq spectrum are related to the occurrence 
of dynamical g-phase transitions and corroborate previ- 
ously obtained numerical results. In Section 6 we express 



ft in terms of the discontinuities of the trajectory scaling 
function a that measures local contraction rates within 
the multifractal attractor. From this relationship we ob- 
tain results for trajectories that involve other multifractal 
regions and show that the dynamics consists of a hierar- 
chy of q-indexes and g-phase transitions. In Section 7 we 
summarize our results and further disscuss the links be- 
tween the thermodynamic and g-statistical schemes and 
call attention to the usefulness of the latter. 



II. THE GOLDEN-MEAN ONSET OF CHAOS 

To assist our presentation in the followingsections we 
recall here basic features of the circle map Q. The 
map is given by 

fn^Kid) =e + n- K/27rsm2TT9, mod 1, (1) 

where the control parameters f2 and K are the bare wind- 
ing number and the amount of nonlinearity, respectively. 
An important parameter for the description of the struc- 
ture of the trajectories of this map is the dressed winding 
number 

(2) 

t^oc t 

the average increment of 6{t) per iteration. The map is 
monotonic and invertible for K < I, develops a cubic 
inflexion point at 6* (or 9—1, see Figs. 1) for K — 1, 
and becomes nonmonotonic and noninvertible for K > 
1. For K < 1 trajectories are periodic (locked motion) 
when w is rational and quasiperiodic (unlocked motion) 
when w is irrational. The winding number u;(f2) forms a 
'devil's staircase' making a step at each rational p/g and 
remaining constant w = p/q for a range of f2. For K — I, 
the critical circle map, locked motion covers the entire f2 
interval leaving only a multifractal set of unlocked fi. The 
threshold to chaos is obtained at the accumulation points 
of sequences of control parameter values of families of 
periodic trajectories. At these accumulation points the 
trajectories are quasiperiodic. For K > 1 the regions 
of periodic motion (Arnold tongues) overlap leading to 
chaotic motion. 

A standard way to study the quasiperiodic route to 
chaos is to fix the value of K = 1 and select an irrational 
w. Then a choice is made of successive rational values 
for to obtain an infinite sequence of winding numbers 
that approach w. A convenient sequence of such ratio- 
nal values for Q are given by consecutive truncations of 
the continued fraction expansion of w. A well-studied, 
an perhaps the most interesting, case is the sequence of 
rational approximants to Wgm = {Vb — l)/2 ~ 0.618034, 
the reciprocal of the golden mean, that yields the wind- 
ing numbers Wn — Fn-i/Fn, where f„ are the Fibonacci 
numbers i^„+i — Fn + F^-i- Within the range of Vl for 
which w = Fn-i/Fn one observes trajectories of period 
Fn, and therefore a route to chaos consists of an infinite 



3 




FIG. 1: Critical circle map K — 1 and supercycle trajectories. 
a.) Q = fl'ao- h) fl = Qoo- The labels F„ correspond to 
trajectory times given by Fibonacci numbers. 



family of periodic orbits with increasing periods of values 
Fn, n ^ oo. It has proved useful to single out from these 
families the periodic trajectories with the 'superstable' 
property. We recall that a superstable orbit 0, 01 is 
an orbit of period T such that df^ {9q) / dO = 0. In our 
case a superstable orbit is one that contains as one of its 
positions 6* = 0. Since the derivative of /a„,i(fi') ai 9 ~ Q 
vanishes the Lyapunov coefficient Ai diverges to — oo for 
these family of orbits. The values of fi, r2„, n = 1, 2, at 
which they occur can be determined from their definition 
when written as = ^n, n = 1, 2, ...The solutions 

of this equation lie on zigzag paths approaching the limit 
f2(x> — 0.606661 which corresponds to Wgm- As the golden 
mean Wgm arises from the condition w = l — w^ , a. second 
superstable family of trajectories is determined by con- 
sidering the winding numbers w'^ = Fn-2/Fn, with Q, 
n = 1,2,..., obtained in turn from /f2/"Y (0) — Fn-i 

converge to VL'^ = 1 — Vloo — 0.393339 which corre- 
sponds to ivjj^ ~ 0.381 966. In Figs, la and lb we show 
the critical map for U, — and U, — f2oo, respectively, 
together with a portion of the superstable trajectory at 
the chaos threshold. In Fig. 2 we show the values for 
the Lyapunov coefficients Ai for if = 1 as a function of 
Q. where the locations of the two superstable families are 
clearly seen and those for their accumulation points are 
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FIG. 2: The Lyapunov coefBcient Ai for the critical circle map 
isT = 1 as a function of winding number. The supercycle di- 
vergences to —00 can be clearly appreciated and the locations 
for the two accumulation points at fi^ and fioo are indicated 
by the arrows. 



indicated. 

A celebrated development Q , |^ dating from the early 
80's is the recognition that the quasiperiodic route to 
chaos displays universal scaling properties. The analysis 
was based on an RG approach analogous to that for the 
period doubling cascade, and can be suitably discussed 
in terms of the critical circle map. The universal prop- 
erties are obtained from the fixed-point map g{9) of an 
RG transformation that consists of functional composi- 
tion and rescaling appropriate for maps with a zero-slope 
cubic inflection point. The fixed-point map satisfies 



g{0) = ag^g[agyng{0/al^), 



(3) 



where agm — —1.288575 is a universal constant 0], 0]. 
This constant describes the scaling of the distance d„ 
from 9 — Q io the nearest element of the orbit with 
e. dn — /n"i(0) — i^n-i. The scaling is given 



by \\mn^ao{.dn I dn+i) = Oigm- Related to the chaos 
threshold at the golden mean Wgm there is a multi- 
fractal set (the attractor) that measures the different 
degrees of clustering of iterate positions. This multi- 
fractal has a dimension spectrum /(a) within a range 
of values between amin and Omax |27l |. The most rar- 
efied region of the multifractal corresponds to amax — 



In? 



/ln|« 



gm\ 



1.897995, and this region is mapped 



onto the multifractal's most crowded region character- 



ized by 



- lnu)g„i/31n \ag 



0.632665 \M- The 



chaos threshold for especial values of the winding number 
is known as local universality (Ij. 
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III. STRUCTURE OF TRAJECTORIES WITHIN 
THE ATTRACTOR 

Wc consider the last of the supcrstablc trajectories, 
with infinite period when fl = floo or fl'^ and start- 
ing at ^(0) = 0. As we shall see all other trajectories 
within the attractor can be obtained from it. We find it 
convenient to plot the positions 9{t) and their iteration 
times t in logarithmic scales as shown in Figs. 3a and 3b, 
obtained for Qoo and fi'^^, respectively, and where the la- 
bels record the times at which positions are reached. A 
distinctive feature in these figures is that positions fall 
along straight diagonal lines, an indication of embedded 
power law behavior. Notice that the positions of the 
main diagonal in Fig. 3a correspond to the times f2n, 
n = 1,2,3,... The succeeding diagonals above it appear 
grouped together. Above the main diagonal there are 
k = 2 close-by diagonals the first with positions that cor- 
respond to times of the form 2F2„ and the second for 
times of the form + -p2n-2- Above these two diago- 
nals there is a group of = 3 diagonals and above them 
a group of A: = 4, etc. The times for the positions along 
each diagonal in the group of k diagonals are, from the 
bottom to the top, of the form kF2n, {k — 1)F2„ + F2n-2, 

{k-2)F2n + 2F2n-2,-, i^2n + (A: - 1)F2„_2 , ^ = 1 , 2, 3, . . . 

The same structure is observed in Fig. 3b with F2n re- 
placed by i^2n-i-i- The slope shared by all the diagonal 
straight lines in Figs. 3a and 3b is — ln|agm| /InWgm. 
Note too that the positions for F2n in Fig. 3a appear 
along the top of Fig. 3b and approach asymptotically 
the horizontal line 6 = 1. And vice versa, the positions 
for ^2n+i in Fig. 3b appear along the top of Fig. 3a 
and approach asymptotically the horizontal line 9 — 1. 
Analogous features hold for the positions of the groups 
of subsequences with A; > 1. A similar pattern of time 
subsequence positions arranged along parallel decreasing 
straight lines is observed in the plot of InjOcc — vs 
Int (or identically. In |f2^ — 9{t)\ vs Int) that is shown in 
Fig. 4. As it can be observed there, the power-law di- 
agonals along which positions fall appear again arranged 
in groups, each with an increasing number of members 
fc = 1, 2, 3, ... The times for the positions along the group 
of k diagonals are, from the bottom to the top, of the form 

kFr, + l,{k- l)^n + K-2 + 1, (fc - 2)Fr, + 2F„_2 + 1,..., 

Fn + {k— l)Fn-2 -I- 1, n = 1, 2, 3, and the slope shared 
by all the diagonals is —3 In \ agm \ / Imugm- 

More specifically, the trajectory starting at 6{Q) = 
maps out the attractor in such a way that the values of 
succeeding (time-shifted r = t + 1) positions 9t form 
subsequences with a common power-law decay of the 
form T-'^/^-Q with Q = l- Inwgm/ In \agm\ ^ 2.897995. 
(To avoid confusion we denote as 6{t) position at time t 
and as 6t the same position at shifted time r). That 
is, the entire attractor can be decomposed into posi- 
tion subsequences generated by the time subsequences 
T = {k — l)Fn + lFn-2, each obtained by running over 
n = 1,2,3,... for fixed values of fc = 1,2,3,... and 
I = 0,1,2, ...,k — 1. Noticeably, the positions for these 
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FIG. 3: a) Positions 0{t) vs t in logarithmic scales for the 
orbit with initial condition 6{0) = at Q'^ of the critical 
circle map K = 1. The lines arc guides to the eye. The 
labels indicate iteration time t, the dashed line and symbols 
X correspond to times of the form 2F2„, and the dotted line 
and symbols -|- correspond to times of the form F2n + -F2n-2- 
The symbols A correspond to the main diagonal positions at 
F2n+i in (b). b) Same as a) with Q'^ replaced by fioo- The 
dashed line and symbols x correspond to times of the form 
2F2n+i, and the dotted line and symbols -|- correspond to 
times of the form F2n+i + -p2n-i- The symbols A correspond 
to the main diagonal positions at F2n in (a). See text for 
description. 



subsequences can be asymptotically obtained from those 
belonging to the superstable periodic orbits of lengths 
Fn at fin or OJj. In particular, the positions for the main 
subsequence k = 1, that constitutes the lower bound of 
the entire trajectory, can be identified to be ~ \dn\ = 
|ctgm|~") f^ = ^n- The positions for each time subse- 
quence that leads a group, fc > 1 and I = 0, i.e. r = kFn 
in Figs. 3 can be expressed, when use is made of the time 
shift T = t + k and 

^"-^ ^^w-"" U) 
F 1 F o"'"Fi ~ ^ ' 

as 

^,=^fcexpQ^(Ag=|t), (5) 

with (5+ = 1 — InWgn, / ln\agm\ and = 
In I agm I /fc In Wgm ■ The same time subsequence positions 
given by r = kFn -|- 1 in Fig. 4 can be expressed as 

\n^-9r\ = \n^-9k\e^PQ_{A^Qh), (6) 
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FIG. 4: Positions \^loc — vs t in logarithmic scales for the 
orbit with initial condition 6(0) = at Qoc of the critical circle 
map K = 1. The labels indicate iteration time t. The dashed 
line and symbols x correspond to times of the form 2F„ + 1, 
and the dotted line and symbols + correspond to times of the 
form Fn + Fn-2 + 1- The same results are obtained with Sloo 
replaced by SI'^q. See text for description. 



with Q_ = 1 — lnwg„i/31n lofgnil and Ag^ — 
'ihi\agm \ /khiWgm- Interestingly, these results can be 
seen to satisfy the dynamical scaling relation, 



h-r = T 



'i/(i"Q) 



0) 



(7) 



with hr = Or or \Q,oo — Or\, where r = kWg^ with Wgm 

3(1-Q_) 
*-gm 1 



" ' , and where g is the fixed- 



point map with = Ok or jfioo — 6k\- The similarities 
with the dynamical properties of the supercycle orbit 2°° 
at the Feigenbaum attractor are remarkable, see Refs. 
0, 0] for the quadratic map and Refs. C3 fo'^ 

generalization to general nonlinearity z > 1. 



IV. SENSITIVITY TO INITIAL CONDITIONS 

We have now elements to derive the sensitivity to ini- 
tial conditions 6(6'o), defined as 6(^(0)) = dO{t)/dO{0), 
when either 0{0) ~ or 0{0) ~ fioo- For simplicity of 
presentation we consider the case of the leading posi- 
tion subsequences within each group k, those for times 
T = kFn as described in the previous section. The results 
for the other subsequences within each group are similar 
and we only quote the differences in the expressions. We 
proceed in the following way. Recall that 



UO) c g{alJ), 



0, 



(8) 



where /„(6') = f^" ^{0) — and consider the lower- 

order terms in the expansion of g^ag^^O), i.e. 



^ 5(0) + g5"'(0)a 



(9) 



If we denote the starting positions of two nearby trajec- 
tories at time r = fc by 



%^h{0)^giO) + -g'"{0)0^ 



and 



0fc=/i(0)~5(O) + i<?"'(O)03, 

their positions at time r — kFn are given by 

1 



and 



ykF„ = UO) c agMO) + -g"\0)af:^0^ 



bkF,, = U<f) ^ a;;j5(0) -I- -g"'{Q)all,f. 



From the above we obtain 

\GkF^ — (j^kF^ I I |2n 



'■gm I 



71 > 1, 



(10) 



(11) 



(12) 



(13) 



(14) 



where we have used the equality sign in consideration of 
the scaling limit for large n. For each subsequence k, the 
sensitivity ^t, given for the trajectories considered by 



UOk) = , Im 



(15) 



can be written, with use of the shifted time variable t 
T — k, and observing that 



'"=fl + ^ 



as the g-exponential 

^4(6*^) = exp^ 



where 



and 



In 



21n| 



0.051003 



t^gm I 



;^(fe) ^ o '^T^\agrn\ 
kin Warn 



(16) 



(17) 



(18) 



(19) 



The value of q in Eq. H18() agrees with that obtained in 
earlier studies (^-jH. 

The departing positions, 0{0) ~ (or 0{0) ~ 1), for 
£_t in Eq. ((T7|l are located in the most crowded region 

of the attractor and the value of k in A^'^' in Eq. H17|) 
locates this position as the top position with time label 
T = kFi — k oi the subsequences described by the times 
T = kFn in Figs. 3a and 3b. By construction, the ob- 
servation times in Eq. (|17|l are precisely r = kFn when 
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the iterate is close to the most sparse region of the at- 
tractor, 9o ~ fioo (or — ^'oo)- The sensitivity ^t+i for 
the inverse trajectories, with departing positions at the 
most sparse region and with observation times when at 
the most crowded region, can be evaluated in the same 
way as above, with the result 



t+i 



exp2-g 



with q as above and where 



(fc) 21n|a 
2-9- kw 



gm 111 '^gm 



(20) 



(21) 



Trajectories expand in one direction and contract in the 
opposite one. We recall that exp (x) = 1/ exp2_ (— x). 



The factor of in 



(k) 



appears because of a basic 
difference between the orbits of periods Fn and Foo- In 
the latter case, to reach 9{t') ~ from 0(0) ~ 1 at times 
t' = Fn -\- \ the iterate necessarily moves into positions 
of the next period Fn+i, and orbit contraction is W^^i 
as effective than expansion. When considering the other 
sets of positions within each group of subsequences one 
obtains the same results as in Eqs. H17|) to Eq. H21|) with 
k replaced by (fc - Z) + in Eqs. I^Hl and (|2J- 



V. MORI'S g-PHASE TRANSITIONS AND 
TSALLIS' q INDEX 

About 15 years ago Mori and coworkers developed a 
broad thermodynamic formalism to characterize sharp 
changes at bifurcations and at other singular attractors in 
low dimensional maps p3 |. The formalism was specially 
adapted to study critical attractors and was applied to 
the specific case of the onset of chaos in the circle map 
at the golden mean 12^, |2^ and in the logistic map at 
the period doubling acctmiulation point [23|, [1^. For 
critical attactors the scheme involves the evaluation of 
fluctuations of the generalized finite-time Lyapunov co- 
efficient 



t-i 
-^Vln 



t > 1, 



(22) 



where f{x) is a map with control parameter(s) tuned 
at a critical attractor. Notice the replacement of the 
customary i by In f in Eq. (|22|l . as the ordinary Lyapunov 
coefficient Ai vanishes for critical attractors at t — > oo. 

The density distribution for the values of A, at i ^ 1, 
P(A,i), is written in the form pl|. [2ll|. ^ 



as 0(q) = — hmt^ooln^(t,q)/lni, where Z{t,(\) is the 
dynamic partition function 



Z[t,q) = dX P{X,t) t-(q-i)^. 



(24) 



The 'coarse-grained' function of generalized Lyapunov 
coefficients A(q) is given by A(q) = d0(q)/dq_and the 
variance v (q) of P{X,t) by t;(q) = dX{q)/dq [H, 
|23 |. The functions (/)(q) and '0(A) are the dynamic coun- 
terparts of the Renyi dimensions D{q) and the spectrum 
/(a) that characterize the geometric structure of the at- 
tractor The form t~ 

4'W for P{X,t) and the weight 
^-(q-i)-'' in Eq. H24|l are intended for the description of 
the envelope of the fluctuations at the critical attactor 

m, p, m. 

As with ordinary thermal 1st order phase transitions, 
a 'g-phase' transition is indicated by a section of lin- 
ear slope = 1 — q in the spectrum (free energy) 
'0(A), a discontinuity at q = q in the Lyapunov func- 
tion (order parameter) A(q), and a divergence at q in 
the variance (susceptibility) v{q). For the onset of chaos 
at golden mean a single (/-phase transition was numeri- 
cally determined |23j| and found to occur approximately 
at a value around rric = — (1 — g) ~ —0.95. It was 
pointed out in Ref. 23] that this value would actually be 
nic = -{1 - q) ^ -amax/2 ~ -0.948997. Our analysis 
below shows that this initial result gives a rough picture 
of the dynamics at the golden mean attractor and that 
actually an infinite family of g-phase transitions of de- 
creasing magnitude takes place there. 



From the results for A^'^^ and Ajj'^ we can construct 
the two-step Lyapunov function 



(fe) 



A, 



(1) 



—oo < q < q, 



A(q) = <0, q<q<2-q, 



(25) 



2 — <? < q < oo. 



where xlj''^ 



21n|Q;gm| /Inwg 



-1.704994. 



1.053744 and A 

5a, 



(1) 



2-9 - 

Direct contact 



See Fig 

can be established now with the formalism developed by 
Mori and coworkers and the g-phase transition reported 
in Ref. |i21||. The step function for A(q) can be integrated 
to obtain the spectrum (j){q) (A(q) = d(j)/dq) and from 
this its Legendre transform ip{X) (= — (1 — q)A). The 
free energy functions 4>{q) and '0(A) obtained from the 
two-step A(q) determined above are given by 



0(q) = < 0, 



A^'i,(q-2 



g < q < 2-g, 

?), q > 2 - q. 



(26) 



and 



p(A,t) =t-'^(^)F(0,t), 



(23) 



where 0'(A) is a concave spectrum of the fluctuations of 
A with minimum 0(0) = and is obtained as the Legen- 
dre transform of the 'free energy' function 0(q), defined 



V'(A) 



q)X, X^i\ < A < 0, 



(1 - q)X, Q<X<X\ 



(1) 



(27) 



See Fig. 5b. The constant slopes of ipiX) represent the q- 
phase transitions associated with trajectories linking two 
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1 .0537 



FIG. 5: g-phase transitions occurring at index values q and 
2 — q. (a) The two-step Lyapunov function A(q). (b) The 
piece- wise spectrum tp{X). These transitions correspond to 
the main discontinuity in the trajectory scaling function o"(y). 
See text for details. 



regions of the attractor, 6* ~ (or ~ 1) and ~ ^oo, 
and their values 1 — q and q—1 correspond to the index q 
obtained for the g-exponentials in Eqs. (|17|l and (|20|l . 
The slope g — 1 ~ — amax/2 — —0.949 coincides with that 
originally detected in Ref. pH. 



VI. MULTIFRACTAL DYNAMICS FROM THE 
TRAJECTORY SCALING FUNCTION 

The trajectory scaling function a{y) for the golden 
mean route is constructed [2^ via the evaluation 

of the ratios 



(T„(m) 



dn 



(28) 



y = Wgrn, 

of course w 



etc. 



and 



grn r u/^^ , lijgm 

gm -rw'g^ = 1. Wc Icave a detailed account of 
these properties to a later occasion . In most cases it 
is only necessary to consider the first few as their rnag- 
nitude decreases rapidly. See, e.g. Fig. 2 in Ref. [2g. 
The discontinuities of cr„(m) can be suitably obtained 
by first generating the superstable orbit at D,oo and then 
plotting the position differences \9r ~ Om \ for times of the 
form T = Fn + m, n — Q, 1,2, in logarithmic scales. 
The distances that separate positions along the time sub- 
sequence correspond to the logarithm of the distances 
dn,m- See Figs. 3 and 4 where the constant spacing of 
positions along the main diagonal provide the values for 
\n.dn,Q — —nhi\agm\ and h\dn,i — — Snlnjagml, respec- 
tively. The constant slope Sm of the resulting time subse- 
quence data is related to an{m), i.e. (t„(to) = {wgmY"'- 
See Figs. 3 and 4 where the slopes of the main diago- 
nal subsequences have the values — In | agm \ / In Wgm and 
— 31n|ag„i| /hiWgm, respectively. From these two slopes 
the value of the largest jump discontinuity of (t„(to) is 
determined. 

The sensitivity ^t(6'(0)) can be evaluated for trajecto- 
ries within the attractor via consideration of the disconti- 
nuities of a{y). We fix the initial and the final separation 
of the trajectories to be the diameters A9{0) = dn.m and 



A0(t) — dn,m+t, 

is obtained as 



t^Fn- 1, respectively. Then, ^((^(O)) 



6(^(0)) = lim 



dn. 



(29) 



Notice that in this limit A6'(0) ^ 0, t — > cxd and the F„- 
supercycle becomes the onset of chaos. Then, ^t(6'(0)) 
can be written as 



CT„(m - 1) 



(t„(to) 



< = F„ — 1 and n large, (30) 



where we have used |cr„(TO)|" ~ n"=i Mi+i,m/rfi,m| and 
di+i^rn+Fn = di+i^rn- Noticc that for the inverse pro- 
cess, starting at A0(O) = dn,m+t — <^n,m-i and ending at 
A9{t') = dn,m = dn,m-i+t', with t' ~ Fn 1 One obtains 



6'(m-l) 



(Tnirn) 



lim - 1) 



<' = F„ + 1 and n large. 



where d„ 



(31) 

Therefore, from the properties of a{y) we can deter- 
mine those for ^j. This function can be obtained at dif- 
ferent levels of approximation by considering only a fixed 
number of its larger discontinuities. The discontinuities 
of cr(2/) appear divided into two sets, the first falls on the 

< 2/ < 1- 



= 0, 1,2..., and where interval < y < Wgm, and the second on 



the denominator is replaced by dn^m — \dm. — dm-F„\ hi 
case m > Fn. The function a{y) is obtained as the limit 
a{y) = lim„^oo cr„(TO) where y = lim„^oom/F„, and 
shows finite (jump) discontinuities at all rationals of the 
form m/Fn- The main discontinuity appears at y = 0, 



The amplitudes of the jumps in the second set are much 
smaller than those in the first set [53 . The simplest ap- 
proximation for cr(y) is to consider only one jump in the 
first set at and none in the second set, so that 



since a{0) = \ag 



but cr(0+) 



and again 



at y 



w 



gm- 



Other discontinuities in a{y) appear at 



< y < w 



2 

gm 7 



<y < 1 



(32) 
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where ao = \agm\ — 2.139583 and ai = \agm\ — 
1.288575. This level of approximation assumes that there 
are only two scaling factors, those corresponding to the 
most crowded and most sparse regions of the attractor. 
The next level of approximation for cr(?/) involves two 
jumps in the first set at 
the second set, so that 



and u;g„j, and still none in 



\ < y < 

(^2^ "^gni <y < Wgm 
"3 \ Wgr,^ <y<l 



y3 



(33) 



where ao = 



ai 



l/a{w 



~ 2.452, a2 = 
l/a{'Wgm) — 1.430 and = \agm\ i22j. The two addi- 
tional scaling factors correspond to two 'midway' regions 
between the most crowded and most sparse regions of 
the attractor. More accurate approximations to a{y) can 
be constructed by incorporation of additional discontinu- 
ities. Note that the number of discontinuities increase in 
each set as Fm — 1, M = 1, 2, ... 

Use of Eq. ^ in Eqs. (EOJl and ^ (together with 
Eq. IjKil) ) recovers our previous results for and Xq in 
Eqs. 1)171) to H21() . and therefore also those for A(q), 0(q) 
and tp{X) in Eqs. to |2ZI). When use is made of 

Eq. H33|) in Eqs. (|30|l and H31|l we obtain three values for 
the q index, qq, qi and q2 (together with the conjugate 
values 2 — go I 2 — and 2 — q2 for the inverse trajec- 
tories). For each value of q there is a set of g-Lyapunov 
coefficients running from a maximum Xq- .max to zero (or 

a minimum A2-qj,min to zero). From the results for A^^ 

(k) 

and X2_q.j j — 0, 1, 2, we can construct three Lyapunov 
functions Aj(q), ~oo < q < oo, each with two jumps lo- 
cated at q = gj = 1 — hi2/ Inaj /aj^i and q = 2 — qj. 
We obtain a couple of two g-phase transitions for each 
of the three values of the q index, go, gi and g2- The 
constant slope values for the g-phase transitions at 1 — go 
and go — 1 appear again, but now we have two other 
pairs of transitions with slope values 1 — gi and gi — 1, 
and, 1 — 92 and 92 — 1, that correspond, respectively, to 
orbits that link the most crowded region of the attrac- 
tor to the 'medium crowded' region, and to orbits that 
link the 'medium sparse' region with the 'most sparse' 
region of the attractor. Similar but more elaborated re- 
sults are obtained when more discontinuities in cr(y) are 
taken into account. These results reflect the multi- region 
nature of the multifractal attractor and the memory re- 
tention of these regions in the dynamics. Thus, when 
a{y) = cr{y+) one has = 1 (or Xq{6{0)) — 0) which 
corresponds to trajectories that depart and arrive in the 
same region, when a{y) ^ <j{y+) the power laws in Eqs. 

and correspond to a departing position in one 
region and arrival at a different region and vice versa, 
the trajectories expand in one sense and contract in the 
other. 



VII. SUMMARY AND ASSESSMENT OF 
APPLICABILITY OF g-STATISTICS 

We have examined in detail the structure of the su- 
perstable orbit at the golden mean critical attractor for 
zero-slope cubic inflection point maps. The knowledge of 
this structure allowed us to determine the sensitivity to 
initial conditions for sets of starting positions within the 
attractor. We found that is made up of a hierarchy of 
families of infinitely many interconnected g-exponentials. 
Each pair of regions in the multifractal attractor, that 
contain the starting and finishing positions of a set of 
trajectories, leads to a family of g-exponentials with a 
fixed value of the index g and an associated spectrum of 

(k) 

g- Lyapunov coefhcients Xq . The indexes g and the spec- 



tra A 



{k) 



and and aI*^' 



■\q come in pairs, q and 2 — g, auu /\q auu 
as these correspond to switching starting and finishing 
trajectory positions. This dynamical organization is diffi- 
cult to resolve from the consideration of a straightforward 
time evolution, i.e. the record of positions at every time 
t for a trajectory started at an arbitrary position 0(0) 
within the attractor. In this case what is observed |2l| 
are strongly fluctuating quantities that persist in time 
with a scrambled pattern structure that exhibits mem- 
ory retention. Unsystematic averages over 0(0) would 
rub out the details of the multiscale dynamical proper- 
ties we uncovered. On the other hand, if specific initial 
positions with known location within the multifractal are 
chosen, and subsequent positions are observed only at 
pre-selected times, when the trajectories visit another se- 
lected region, a distinct g-exponential expression for is 
obtained. 

Subsequently, we explained that the above described 
properties for the sensitivity relate to sets of dynamical 
phase transitions. That is, the dynamics on the attractor 
consists of an infinite family of Mori's g-phase transitions, 
each associated, as before, to trajectories that have com- 
mon starting and finishing positions located at specific 
regions of the attractor. The specific values of the vari- 
able q in the formalism of Mori and colleagues at which 
the g-phase transitions take place are the same values for 
the Tsallis entropic index g in ^4. The transitions come 
in pairs at g and 2 — g as they are tied down to the expres- 

ik\ 

sions for Xq in ^j. The dominant dynamical transition 
is associated with movement from the most crowded to 
the most sparse regions of the attractor (or vice versa) 
and the values of its relevant parameters coincide with 
those found in an earlier study with the use of Mori's 
formalism 21]. Here we have determined the g-phase 

transitions via the direct evaluation of the spectra A^*^^ 

and }y2-q Also, we established a connection between 

the trajectory scaling function u and the expressions for 
^t, such that the discontinuities of the former determine 
the g-phase transitions and the values for the g index. 
Conversely, from the g-exponential expressions for the 
scaling function a can be conveniently determined. 
Our results clearly apply to many other families of 
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maps with zero slope inflection points of cubic nonlin- 
earity and are likely to hold also for general nonlinearity 
z > 1 These results may be of relevance in con- 

densed matter problems where the quasiperiodic route 
to chaos is involved. One interesting example is the phe- 
nomenon of the localization transition for transport in 
incommensurate systems, where Schrodinger equations 
with quasiperiodic potentials [sj are equivalent to non- 
linear maps with a quasiperiodic route to chaos, and 
where the divergence of the localization length translates 
into the vanishing of the ordinary Lyapunov coefficient 
[3^ . The dynamics at this type of multifractal critical 
attractor is also represented by a given by families of 
interconnected g-exponentials with g-generalized spectra 

(k) 

of Lyapunov coefficients Xq . The latter quantities yield 
information on the power-law decay rates of multifractal- 
structured wave functions at the localization transition. 
Our flndings reported here are very similar to those re- 
cently determined for the Feigenbaum attractor in uni- 
modal maps and reinforce the notion of generality 
of our identification of the source for the entropic index 
q observed at critical attractors. 



subsequence T{n,k,l), as explained in Sections 3 and 4 
and shown in Figs. 3 and 4. For simplicity we consider 
now only the case I — 0, i.e. the subsequences r(n, k) = 
kFn, n = 1,2,3,... for fixed values of fc = 1,2,3,..., but 
it is important to note that every natural number r ap- 
pears in the complete set of subsequences. One branch 
of the envelope of (_t with 6*0 = is visited at times t 
of the form r = F„ for which we find A(F„,0) = A^^-* 
while for other time subsequences one has A(A:_F'„,0) = 
Ag'^'' (see Eq. (I19II . Again for simplicity we ignore the 



(see 



inverse trajectories that lead to the conjugate A2_g 
Eq. (jai). 

As we have proved, along each time subsequence 
the sensitivity — \df''^\0o)/d9o\ is given by a q- 
exponential, like in Eq. IjlSI) . and consequently we can 
express X{t,9o) in terms of the g-deformed logarithm 
A(t, Oo) = In, |d/(")(0o)/rf^o| with q given by e.g. Eq. 
(|18|l . Therefore for times of the form r(n, k) — kFn we 
have 

P(A, r)Wi\, t) = 5{X - Af )) exp,(A«i), (36) 
with t = T — k, and 



A. The utility of q-statistics 

Before concluding our analysis it is of interest to recon- 
sider our developments in such a way as to illustrate the 
incidence and usefulness of the q-statistics in the descrip- 
tion of the dynamics on the attractor at the transition to 
chaos. We return to the dynamical partition function in 
Eq. H24(l and rewrite it as 



Z(T,q)= / dAP(A,T) WiX,T) 



i-q 



(34) 



where the weight W{X, r) was assumed by Mori and col- 
leagues [i^ to have the form W{X,t) = exp(AT) for 
chaotic attractors and W{X, r) = for the envelope of 
fluctuations of the attractor at the onset of chaos, with A 
a generalized Lyapunov coefficient. The variable q plays 
the role of a 'thermodynamic field' like the external mag- 
netic field in a thermal magnet or, alternatively, 1 — q can 
be thought analogous to the inverse temperature. We 
note that W{X, r) can be written in both cases as 



W{X,t) 



d0o 



(35) 



since X{T,eo) ^ A{t) In \df^^H9o)/deo\ with A{t) = r^i 
for chaotic attractors and A{t) — (Inr)^^ for the enve- 
lope of the fluctuating sensitivity of the critical attractor 
[23 ■ For chaotic attractors the fluctuations of A(r, 0o) 
die out as r ^ 00 and this quantity becomes the ordi- 
nary Lyapunov coefficient Aoo independent of the initial 
position 60 . For critical attactors A(t, ^o) maintains its 
dependence on r and 9q but as we have seen A(r, Sq) is a 
well-defined constant provided r takes values from a time 



Z(t,q) = iy(AW,i)i- 



(l-q)/(l-q) 



(37) 

Notice that in Eq. 1)37(1 q is a running variable while q is 
fixed as in Eq. fll^ . 

The next step is to consider the uniform probability 



distribution, for fixed Xq'^^ and t, given by pi{Xq'^' ,t) = 

W ,1 = 1,..., W, where W is the integer nearest to (a 
large) W. A trajectory starting from 9o = visits a new 
site 9r belonging to the time subsequence r = fcF„ every 
time the variable n increases by one unit and never re- 
peats one. Each time this event occurs phase space is pro- 
gressively covered with an interval A9n = 9kF„+i ~ ^^kF^ ■ 
The difference of the logarithms of the times between 
any two such consecutive events is the constant Inw"^, 
n large, whereas difference of the logarithms of the cor- 
responding phase space distances 9kF„+i and 9kF„ is the 
constant 21n|agm|- We define Wn to be the total phase 
space distance covered by these events up to time r = 
kFn- In logarithmic scales this is a linear growth pro- 
cess in space and time from which the constant prob- 
ability Wn of the uniform distribution pi is defined. 
In ordinary phase space 9 and time t we obtain instead 
= exp(-Af ^t). We then have that 

w 

z(t,q) = ^pHi + (i-q)^q, 

i=l 



Ak) 



(38) 



where S'q = hiq W is the Tsallis entropy for pi. 

The usefulness of the g-statistical approach is now ev- 
ident when we recall from our discussion in the previ- 
ous section that the dynamics on the critical attractor is 
constituted (in its entirety) by a discrete set of g-phase 
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transitions and that at each of them the field variable 
q takes a specific value q = q. At each such transition 
both Z and Sq grow linearly with time along the subse- 
quences T(n, k). In addition, Ag'"'' can be determined from 
Sq/t — Xq^\ Thus, with the knowledge we have gained. 



a convenient procedure for determining the fluctuating 
dynamics at a critical multifractal attractor could take 
advantage of the above properties. 
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